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1) ■ Abstract 

^ , We develop a non-life reserving model using a stable-1/2 random bridge to 

simulate the accumulation of paid claims, allowing for an arbitrary choice of a 
priori distribution for the ultimate loss. Taking a Bayesian approach to the 
reserving problem, we derive the process of the conditional distribution of the 
ultimate loss. The 'best-estimate ultimate loss process' is given by the condi- 
tional expectation of the ultimate loss. We derive explicit expressions for the 
best-estimate ultimate loss process, and for expected recoveries arising from ag- 
gregate excess-of-loss reinsurance treaties. Use of a deterministic time change 
allows for the matching of any initial (increasing) development pattern for the 
paid claims. We show that these methods are well-suited to the modelling of 
claims where there is a non-trivial probability of catastrophic loss. The gener- 
alized inverse- Gaussian (GIG) distribution is shown to be a natural choice for 
the a priori ultimate loss distribution. For particular GIG parameter choices, 
the best-estimate ultimate loss process can be written as a rational function of 
the paid-claims process. We extend the model to include a second paid-claims 
process, and allow the two processes to be dependent. The results obtained can 
be applied to the modelling of multiple lines of business or multiple origin years. 
The multi-dimensional model has the property that the dimensionality of cal- 
culations remains low, regardless of the number of paid-claims processes. An 

S^ algorithm is provided for the simulation of the paid-claims processes. 
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1 Introduction 

The term 'stable process' refers here to a strictly stable process with index a £ (0, 2); 
thus, we are excluding the case of Brownian motion (a = 2). The use of stable processes 
for the modelling of prices in financial markets was proposed by Mandelbrot [26( in 
connection with his analysis of cotton futures. In the tails, the Levy densities of stable 
processes exhibit power-law decay As a result, the behaviour of stable processes is 
wild, and their trajectories exhibit frequent large jumps. The variance of a stable 
random variable is infinite. If a < 1, the expectation either does not exist or is infinite. 
This heavy-tailed behaviour makes stable processes ill-suited to some applications in 
finance, such as forecasting and option pricing. 

To overcome some of the drawbacks of the stable processes, so-called tempered 



stable processes have been introduced (see Cont & Tankov [14j, for example, for details). 
A tempered stable process is a pure-jump Levy process, and its Levy density is the 
exponentially dampened Levy density of a stable process. The exponential dampening 
of the Levy density improves the integrability of the process to the extent that all the 
moments of a tempered stable process exist. Tempered stable processes do not possess 
the time-scaling property of stable processes. 

In this paper we apply stable- 1/2 random bridges to the modelling of cumulative 
losses. The techniques presented can equally be applied to cumulative gains. The 
integrability of a stable-1/2 random bridge depends on the integrability of its terminal 
distribution. At some fixed future time, the nth moment of the process is finite if 
and only if the nth moment of its terminal value is finite. Thus a stable-1/2 random 
bridge with an integrable terminal distribution can be considered to be a dampened 
stable-1/2 subordinator. In fact, the stable-1/2 random bridge is a generalisation of 
the tempered stable-1/2 subordinator. If the Levy density of a stable-1/2 subordinator 
is exponentially dampened, the resulting process is an IG process. We shall see that 
the IG process is a special case of a stable-1/2 random bridge. 

We look in detail at the non-life reserving problem. An insurance company will incur 
losses when certain events occur. An event may be, for example, a period of high wind, 
the flooding of a river, or a motor accident. The losses are the costs associated with 
recompensing policy-holders who have been disadvantaged by an event. These costs 
might, for example, cover repairs to property, replacement of damaged items, loss of 
business, medical care, and so on. Although a loss is incurred by the insurance company 
on the date of an event (the 'loss date' or 'accident date'), payment is generally not 
made immediately. Delays will occur because loss is not always immediately reported 
to the company, because the full extent of the costs takes time to emerge, because the 
insurance company's obligation to pay takes time to establish, and so forth. 

In return for covering policy-holder risk, the insurance company receives premiums. 
The premiums received over a fixed period should, typically, be sufficient to cover the 
losses the company incurs over that period. Since losses can take years to pay in full, 
the company sets aside some of the premiums to cover future payments; these are 
called 'reserves'. If the reserves are set too low, the company may struggle to cover its 



liabilities, leading to insolvency. Large upward moves in the reserves due to a worsening 
in the expected future development of liabilities can cause similar problems. If the 
reserves are set too high, the company may be accused by shareholders or regulators of 
inappropriately withholding profits. Hence, it is the interest of the company to forecast 
its ultimate liability as accurately as possible when deciding the level of reserves to set. 

We use a stable-1/2 random bridge to model the paid-claims process (i.e. cumu- 
lative amount paid to-date) of an insurance company. The losses contributing to the 
paid-claims process are assumed to have occurred in a fixed interval of time. Sometimes 
claims-handling information about individual losses is known, such as that contained 
in police or loss-adjuster reports. In the model presented here, such information is 
disregarded, and the paid-claims process is regarded as providing all relevant infor- 
mation. We derive the conditional distribution of the company's total liability given 
the paid-claims process. We then estimate recoveries from reinsurance treaties on the 
total liability. The expressions arising in such estimates are similar to the expectations 
encountered in the pricing of call spreads on stock prices. 

We shall examine the upper-tail of the conditional distribution of the ultimate 
liability, and find that it is as heavy as the a priori tail. This has an interesting 
interpretation in the case when the insurer is exposed to a catastrophic loss. At time 
t <T, the probability of a catastrophic loss occurring in the interval [t, T] decreases as t 
approaches T. However, in some sense, the size of a catastrophic loss does not decrease 
as t approaches T, since the tail of the conditional distribution of the cumulative loss 
does not thin. 

When the a priori total loss distribution is a generalized inverse-Gaussian distri- 
bution, we find that the model is particularly tractable. We present a family of special 
cases where the expected total loss can be expressed as a rational function of the current 
value of the paid-claims process. That is, each member of the family is a martingale 
that can be written as a rational function of an increasing process. 

The model can be extended to include more than one paid-claims process. We 
consider the case where there are two processes that are not independent, and which 
have different activity parameters. We then have two ultimate losses to estimate. We 
provide expressions for the expected values of the ultimate losses given both paid- 
claims processes. The numerical computations required to evaluate these expressions 
are no more difficult than those of the one-dimensional case. We demonstrate how to 
calculate the a priori correlation between the ultimate liabilities. The correlation can 
be used as a calibration tool when modelling cumulative losses arising from related 
lines of business (e.g. personal motor and commercial motor business). 

We also describe how to simulate sample paths of the stable-1/2 random bridge, 
and how to use a deterministic time-change to adjust the model when the paid-claims 
process is expected to develop non-linearly. 



2 Preliminaries 

We fix a probability space (fl, Q, J 7 ), and assume that all processes and nitrations under 
consideration are cadlag. Unless otherwise stated, when discussing a stochastic process 
we assume that the process takes values in R and begins at time 0, and the filtration 
is that generated by the process itself. We work with a finite time horizon [0, T]. 

2.1 Levy processes and stable processes 

A Levy process is a stochastically-continuous process that starts from the value 0, and 
has stationary, independent increments (see Bertoin [7(], Kyprianou 24|, and Sato 32J). 



An increasing Levy process is called a 'subordinator'. If a Levy process {5" t Q } satisfies 
the scaling property 

{r 1/a £ fet } t >o = {S t } t > for k > 0, (1) 

then we say that it is a (strictly) stable processes with index a, or a stable-a process. 
It can be shown that for ([T]) to hold we must have a G (0,2]. A stable-2 process is a 
scaled Weiner process, and a stable- 1 process is a Cauchy process. 

If {Sf} is a stable subordinator, the Laplace transform of 5° exists and is given by 

E[e" AS * Q ] = exp{-Kt\ a ) for A > 0, (2) 

where k > 0, and a must be further restricted to a e (0, 1). 

2.2 Stable-1/2 subordinator 

Let {St} be a stable-1/2 subordinator . The Laplace transform of St is 

E[e-^]=exp|-^| (A>0), (3) 

for some c > 0. Then {S t } satisfies the scaling property {A; -2 ^} = {St}, for k > 0. 
The random variable St has a 'Levy distribution' with density 

ft{x) = t{x>0} 7T^ exv {-2— l (4) 



We call c the 'activity parameter' of {S t }. The density @ is bounded for all t > and 
is strictly positive for all x > 0. Integrating (|3j) yields the distribution function 



/ f t (y)dy = 2^[-ctx- 1 / 2 } 
Jo 



(5) 



where $[x] is the standard normal distribution function. The random variable St has 
infinite mean, indeed E[5f] < oo if and only if p < 1/2. The density of 1/S t is 



x M- 1; 



cf 



-,/ _1 ' 2 exp 
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(6) 



{a;>0} V2r[i/2] 

Thus the increments of {S t } are distributed as reciprocals of gamma random variables. 
For {W t } a Weiner process, define the exceedence times {r t } t > by 

r t = mi{s : W s > ct}. (7) 

From Feller |l9l X.7], we then have 



law 



{S t } = {n}. 



(8) 



3 Stable-1/2 random bridge 



Levy random bridges are defined in Hoyle et al. [22[ , where some of their properties are 
derived. We call the random bridge of a stable-1/2 subordinator a 'stable-1/2 random 
bridge'. We shall state results about stable-1/2 random bridges, and refer to 22J for 
further details. 



3.1 Stable-1/2 bridge 

Fix z > and let {S^}o<t<T be a bridge of the process {St} to the value z and time 
T. That is, {S^ } is a stable-1/2 subordinator conditioned to arrive (and terminate) 
at z at time T. Levy bridges are Markov processes, and the transition law of {S^ } is 



j(«) 



S^ e dy 
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(y-x;z-x) dy, 



where 



fa(y;z) 



ft{y)fr-t(z - y) 



Hz) 



i 



{0<y<z}- 



1 ct(T-t) ex P 



1 c 2 (Ty-tz) 2 

2 yz(z-y) 



^7T 



T 



(y - y 2 /zf' 2 



Note that ftriv] z ) is the density function of the random variable S 1 
bounded, and has bounded support, so 



tT • 



(9) 

(10) 

(11) 
This density is 



E 






< oo for p > 0. 



(12) 



Integrating the density ffTTl) yields the following distribution function for y G [0, z\: 



F tT {y;z) = <5> 



c{Ty - tz) 
\/yz{z - y) 



2t ^ e 2ch(T-t)/z $ 

T ' 



c{{2t-T)y-tz) 
y/yzjz - y) 



(13) 



Remark 3.1. When t = T/2, the second term in the distribution function [W\l van- 
ishes. The distribution function is analytically invertible, and we obtain the following 
identity, where Z is a standard normal random variable: 



q {z) law 1 , -, 
°T/2,T ~ O ' 



^c 2 T 2 /z + Z 2 



(14) 



Proposition 3.2. For fixed k > 0, the stable-1/2 bridge {S t ^ } satisfies the scaling 
property 

{*®\.„ = {*-*&*_ . (15) 



0<t<T I ' J 0<t<T 

Proof. The transition probabilities of the processes are given by 



S ,(z) < v 



^sT — X 



i-2 o(k 2 z) ^ 



^ °ks,kT ~ x 



F t _ S)T _ s (y- x;z- x), 
F k ( t -s),k{T-s)(k 2 y - k 2 x; k 2 z - k 2 x) 



(16) 
(17) 



for < s < t < T. Substituting for Ftriy] z) given in ([TBI shows that these probabilities 
are equal. 

Proposition 3.3. 



{Stf ■ } — > {^ z } as c —)• oo. 
Proof. Fix z > 0. It is sufficient to show that 

lim F tT (y; z) = t {Ty>tz} for Lebesgue-a.e. y G (0, z) 

C— 7>CO 

since this is equivalent to 

< e 



lim ( 

c— »oo 



Q tT rpZ 



Define a by 



1 for all t G [0, T] and any e > 0. 
(2* - r)y - tz 



\/yz{z - y) 
and note that a > for y G (0, z). The inequality (l|, 7.1.13] states that 



e-* d£ < 



x 



+ v/x 2 + 4/ 



7T 



from which we deduce 



e 2c2 < T -*)/^[-ac] < e 



2c 2 t(T-t)/z 



{X > 0), 



,-a 2 c 2 /2 



+ v/« 2 c 2 + 2/tt 



TTac 



2 exp i c ^(z-y)* 



vrac + v /a 2 c 2 + 2/vr' 
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(18) 

(19) 

(20) 
(21) 
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(23) 



Since the left-hand side of ( 1231) is positive, we see that 

lim e^W-^Ql-ac} = 0. 



(24) 



Then we have 

lim F t T{y; z) = lim $ 



c(Ty - tz) 



!_^j i ime 2A(T-t)/ z$[ _ ac ] 

T J c— »oo 



JUTi,-**>m. — ttIj 



y/yz(z--v) 

{Ty-tz>0} ~ ^-{Ty^tz}, 

which completes the proof. 

We define the incomplete first moment Mtriy', z) of S^ by 

M tT (y, z)= f u f tT (u; z)du (0<y<z). 
Jo 

Straightforward use of calculus gives 



M tT (y; z) = -z{<$> 



c{Ty - tz) 
^yz(z - y) 



^<?t{T-t)/z $ 



?(*) 



c{{2t-T)y-tz) 
^yz(z - y) 



We can also calculate the second moment of S tT . The result is 
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It then follows from equation ([9]) that for < s < t < T we have 
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Figure 1: Simulations of the stable-1/2 bridge demonstrating the influence of the ac- 
tivity parameter c. Qualitatively speaking, increasing the value of c decreases the 
frequency of large jumps, and increases the frequency of small jumps. 

3.2 Stable-1/2 random bridge 

Let v be a probability law on IR + . We say that {£tr}o<t<T is a stable-1/2 random 
bridge with terminal law v if the following conditions are satisfied: 

1. £tt has law v. 

2. There exists a stable-1/2 subordinator {St} such that 

Q [6i,t <zi,..., 6„,t < £ n | £tt = z] = Q [S tl < xi, . . . , Stn < x n | S T = z] , 

(31) 
for every n G N + , every < t 1 < ■ ■ ■ < t n < T, every (xi, . . . ,x n ) G W 1 , and 
u-a.e. 2. 

Hence, if the bridge laws of {£tr} are the bridge laws of a stable-1/2 subordinator, then 
{£ir} is a stable-1/2 random bridge. It is useful to think of {£*t} as a stable-1/2 bridge 
to a random variable with law v. 

The finite-dimensional distributions of {£tr} are 

n 

Q [£ tl , T e dxi, . . . ,6„,t e dx n ,^ TT G dz] = J [ [/ti-ti-i (a?< - «i-i) da*] ^ tn (dz;a; n ), 

»=i 

(32) 

where the (un-normalised) measure ?/>t(dz;£) is given by 
7/>o(dz;0 = u(dz), 



(33) 



/t(^) 



4*>«> 



(i-9 



n3/ 2 ex P| 2 



1 /c 2 T 2 c 2 (T-tf 



z-i 



u(dz), (34) 



for < t < T. It can be shown that {£tr} is a Markov process with transition law 



tr e dy | £ sT = x] 
tt e dy | £ sT = a?] 



ft- s (y-x)dy, 



V> 8 (dy;g) 
^ s (R;z)' 

for < s < t < T. 

Now fix a time s < T and define a process {t]tT}s<t<T by 

Then {??4t} is a stable-1/2 bridge with terminal law fT- s (x)ipT-& 
more, given £ s j<, {^t} is a stable-1/2 bridge with terminal law 



(35) 



v*{A) 






(36) 
x) dx. Further- 

(37) 



where A + y denotes the set {x : x — y G A}. We call this the 'dynamic consistency' 
property of stable-1/2 bridges. The financial significance of the dynamic consistency 
condition is discussed in [9j , [ll| , and 22| . 



4 Insurance model 

We approach the non-life insurance claims reserving problem by modelling a paid-claims 
process by a stable-1/2 random bridge. We shall look at the problem of calculating 
the reserves required to cover the losses arising fro m a sing le line of business when 
we observe the paid-claims process. Arjas |2j and iNorbergl 3(1 l31| provide general 
descriptions of the problem. England & Verrall [18| and Wurthrich & Merz 33j survey 



some of the existing actuarial models. Biihlmann |l2j and Mikosch [29| contain related 



topics. The present work ties in with that of Brody et al. [TJJ who use a gamma random 
bridge process to model a cumulative loss or gain. 

The method we use has a flavour of the Bornhuetter- Ferguson model from actuarial 
science [8] (see also [18j). In implementing the Bornhuetter- Ferguson model, one begins 
with an a priori estimate for the ultimate loss (the total cumulative loss arising from 
the underwritten risks). Periodically, this estimate is revised using a chain-ladder 
technique to take into account the a priori estimate and the development of the total 
paid (or reported) claims to date. 

In the proposed model, we assume an a priori distribution for the ultimate loss. 
By conditioning on the development of the paid-claims process, we revise the ultimate 
loss distribution using Bayesian methods. In this way we continuously update the 
conditional distribution for the total loss. This is as opposed to the deterministic 
Bornhuetter-Ferguson model in which only a point estimate is updated. Knowledge 



of the conditional distribution allows one to calculate confidence intervals around the 
expected loss, and to calculate expected reinsurance recoveries. 

Credibility theory uses Bayesian methods to calculate insurance premiums (see, for 
example, Buhlmann & Gisler [13|). In a typical set-up, the premium a policyholder 
must pay for insurance cover is a functional of their future claims distribution. The 
policyholder's future-claims distribution is parameterised by a rating factor B, whose 
value is unknown, but for which we have an a priori distribution. Observation of 
the policyholder's claims history then leads to an updating of the distribution of 



and hence an updating of their premium for future cover. Wiirthrich & Merz [33 
discuss Bayesian reserving models based on credibility theory. The two main differences 
between credibility reserving models and the present work are (i) they only update in 
discrete time, and (ii) the ultimate-loss distribution is updated indirectly through a 
'rating' variable G. The second difference means that it is not straightforward to allow 
an arbitrary a priori ultimate-loss distribution. 

The main assumptions of the stable- 1/2 bridge model are the following: 

1. The claims arising from the line of business have run off at time T. That is, at 
time T all claims have been settled, and the ultimate loss Ut is known. 

2. Ut has a priori law v such that Ut > and E[C/|.] < oo. 

3. The paid-claims process {£*t} is a stable-1/2 random bridge, and £tt — Ut- 



4. The best estimate of the ultimate loss is U t T = E 
natural filtration of {£tr}- 



Un 



•'t 



where {J-jf } is the 



A few remarks should be made about assumption HI First, using the natural filtration 
of {£tr} as the reserving filtration means that the paid-claims process is the only source 
of information about the ultimate loss once the measure v is set. We do not consider the 
situation where we have access to information about claims that have been reported 
but not yet paid in full (such as case estimates). Second, the expectation is taken 
with respect to Q, which may or may not be the 'real-world' measure. Let us call Q 
the actuarial measure. When reserving, practitioners routinely discount data before 
modelling. Discounting may adjust the data for the time-value of money or for the 
effects of claims inflation. Claims inflation, and interest rates, though understood to be 
stochastic, usually only provide a small amount of uncertainty to the forecasting of the 
ultimate loss, relative to the uncertainty surrounding the frequency and (discounted) 
sizes of insurance claims. Furthermore, it is often for practical purposes reasonable 
to assume that claims inflation and interest rates are independent of claim frequency 
and size. Hence, a stochastic reserving model may lose little from the assumption 
that interest rates and inflation rates are deterministic. We make this assumption, 
and further assume that the paid-claims process has been discounted for the effects of 
interest and inflation. 
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5 Estimating the ultimate loss 



The time-t conditional law of Ut is 



u t (dz) 



M^&t) 



1 



{zHtr} 



3/2 



(T-t) 



ex P -T ^"T H^ 



3/2 



The best-estimate ultimate loss is then 






T 2 



(38) 



£ M*0 



u 



1T 



zv t (dz) 



(39) 



Ot 



At time t < T, the total amount of claims yet to be paid is Ut — £,tr- The amount that 
the insurance company sets aside to cover this unknown amount is called the reserve. 
The expectation of the total future payments is called the best-estimate reserve, and 
can be expressed by 

RtT = U tT - itT- (40) 

For prudence, the reserve may be greater than the best-estimate reserve. However, for 
regulatory reasons it is sometimes required that the best-estimate reserve is reported. 
The variance of the total future payments is the variance of the ultimate loss, which is 



Var 



u T -C, 



tr 



tS 



Var 



Un 



JF} 



[z-U 



tT 



u t (dz). 



(41) 



itr 



6 The paid-claims process 



We give the first two conditional moments of the paid-claims process. Using equations 
and f !30p . and a straightforward conditioning argument, we have 



e[6t|^] 



T-t t-s 

-£sT + 7^ ^sT, 



T 



T 



(42) 



and 



E[&] = % f A 1 -c(T-t)e^V^$ 1-cTz-^] \ v& 



T 



t 
T 



-E [U T ] -c(T-t)V2% I z 3/2 e c — $ [-cTz~ 1/2 } v(dz). (43) 



Equation (1421) implies that the paid-claims development is expected to be linear. We 
return to this point later. Fix s < T and define the relocated process {rj t T}s<t<T by 



VtT — £tr — CsT- 



(44) 
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The dynamic consistency property implies that, given £ s t, {rjtr} is a stable-1/2 random 
bridge with marginal law of tjtt being v*(A) = v s (A + £ s t)- Then we have 



E[^t | Jl] = %'t I &r] + 2&rE [ V tr I U] + l 



2 
sT 



T 



T 



~A + 



c(T - 1) 



T 

2~TT 



-E[U T \£ sT } 



6t) 5 e 2 < z -«-T) $ 



c(T - a) 



v.(dz). (45) 



7 Reinsurance 

An insurance company may buy reinsurance to protect against adverse claim develop- 
ment. The stop-loss and aggregate excess-of-loss treaties are two types of reinsurance 
that cover some or all of the total amount of claims paid over a fixed threshold. Under 
a stop-loss treaty, the reinsurance covers all the losses above a prespecified level. If 
this level is K, then the reinsurance provider pays a total amount (Ut — K) + to the 
insurance company. The 'aggregate L excess of K' treaty is a capped stop-loss, and 
covers the layer [K, K + L] . In this case the reinsurance provider pays an amount 
(U T - K) + -(U T -K- L) + . 

The insurance company typically receives money from the reinsurance provider 
periodically. The amount received depends on the amount that they have paid on 
claims to-date. If the insurer has the paid-claims process {£,tr}, and receives payments 
from a stop-loss treaty (at level K) on the fixed dates t\ < t<i < ■ ■ ■ < t n = T, then the 
amount received on date £,■ is 



K) + . (46) 

16]) can be calculated using the 



(&,t - K)+ - fo_ lfT 

The expected value of reinsurance payments such as 
following: 

Proposition 7.1. Fix t G (0, T). At time s <t, the expected exceedence of £tr over 
some fixed K > is 

D st = E[(fa-K) + \j*] 



T-t t~s TT 

T — s T — s 



K 



/■oo 

+ t { K HsT }(K - £ sT ) / F t . s ^ s (K - i 

J K 

M t _ SjT _ s (K -£ sT ;z- £ sT 



K 



Proof If K < £ sT then 

E [(&r - K)+ | J*] = E [&r | J*] - K 

T-t, t- 

= 6t + 



, T ]z - Z, sT ) v s (dz) 
u g (dz). 



(47) 



T-s 



U sT - K. 



(4£ 
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Thus we need only consider the case when K > £ sT . The J-"|-conditional law of ^ t r is 

tT edy\^} = ^§^\fUy-U)dy. (49) 



Hence we have 



1 f°° 

Dst = / nn> t \ / & ~ ^)^( R ! V)ft-s(y ~ 6r) dy 
ip s (R;£ sT ) J K 

' r (y - K) r fr ^ ~ y) v(dz) ft- s (y - U) dy 



ip 3 (R;£ sT ) J K J K f T (z 



{y _ K) fr-«(* ~ v)Mv ~ M dy u(dz) 



*/j s (R;£ sT ) J K J K f T (z 

/*oo pz 

= / {y-K)f t _ Sj T- s {y-£ S T-,z-£ sT )dyv s {<iz 

JK JK 

Making the change of variable x = y — £ s t yields 

D 8t = I I {x + £sT ~ K)f t -sF- 9 {x', z - i sT ) dx v s {dz) 

Jk JK-e,T 



(50) 



K JK-( sT 
°° f t - S 1 

<^ YZT S ( Z ~ ^t) - M t _ StT _ s (K - £st; z - i sT ) \ v s {dz) 

POO 

+ (U -K) {1 - F t ^ T ^ s (K - £ sT ; z - £ sT )} v,(dz) 

JK 

t-s C°° 

U + ^ U sT -K + {K- £ sT )F t -s,T-s{K - Ul z - U) v s (dz) 

T-s J K 



T-t u t-s 



T-s 

"0O 

M t _ s , T _ s (K - £ sT ; z - £ sT ) is s (dz). (51) 

>K 

D 

Suppose that the insurance company has limited its liability by entering into a 
stop-loss reinsurance contract. At time s G [0, T), the expected reinsurance recovery 
between times t and u is 

E [(U - K) + - (6t - K) + | J*] = D su - D st , (52) 

for s < t < u < T. 

Using a similar method to the calculation of D st , we can calculate the expectation 
of £tT conditional on it exceeding a threshold. For a threshold 9 > £ s t, we find 

m ,, p >fll ^sT + j^UsT ~ C ^j - 6t; * - fej ^(dz) 

Sometimes called the conditional value-at-risk (CVaR), this expected value is a coherent 



risk measure, and is a useful tool for risk management (see McNeil et al. J28j]). Note 
that CVaR is normally defined as an expected value conditional on a shortfall in profit. 
Since we are modelling loss, and not profit, the risk we most wish to manage is on the 
upside. Hence, conditioning on an exceedence is of greater interest. 
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8 Tail behaviour 

In this section we consider how the probability of extreme events is affected by the paid- 
claims development. Suppose that the line of business we are modelling is exposed to 
rare but 'catastrophic' large loss events. In this case we assume that the a priori 
distribution of the ultimate loss has a heavy right-tail. If a catastrophic loss could 
hit the insurance company at any time before run-off, then it is important that any 
conditional distributions for the ultimate loss retain the heavy-tail property. We shall 
see that in the stable-1/2 random bridge model, the conditional distributions are as 
heavy tailed as the a priori distribution. 

Assume that Ut has a continuous density p(z) which is positive for all z above some 
threshold. Then the value of Ut is unbounded in the sense that 



T T > x] > 0, for all igR. (54) 

Define 

Tail, = li,„ n[ , QB 7 >J :' r (55) 

l^oo y [£_ TT - £ tT > L | t,tT\ 

If Taih = oo then the tail of the future-payments distribution at time t > is not as 
heavy as the a priori tail. That is, a catastrophic loss at time t is 'smaller' than a 
catastrophic loss at time 0. If Tail t = then the tail of the future-payments distribution 
is greater at time t than a priori. If < Tail t < oo then the tail is as heavy at time t 
as a priori. Using L'Hopital's rule, we have 



Taih = lim 



«R;&r)jrpM<i2 



^/^U^r^B (££-?))*)<* 



lim 



*(«;««■) p(l) 



= *«^£i^&- < 56) 

for t G (0, T). Some examples follow: 

• If p(z) oc l{ 2>0 }e _z (exponential) then Taih = ipt(M', £tr) <^ tT ■ 

• If p(z) oc l{ z> o}e _z (half-normal) then TaiL, = ipt{^',£tT) ^ tT - 

• lfp(z) oc li^oyz-We- 1 '* (Levy) then Tail* = ^ t (R;&r). 

This property has an interesting parallel with the subexponential distributions. By 
definition, X has a subexponential distribution if 

hm —7-7 — = n, (57) 
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where {Xj}" =1 are independent copies of X (see Embrechts et al. |17j). We note that 



r Q[Zt>L] 

££oQ[Z T -Z t >L\Z t ] =0 °> (58) 

for {Z t } a Brownian motion, a geometric Brownian motion or a gamma process. If 
{Z t } is a stable-1/2 subordinator, so the increments of {Z t } are subexponential, then 

lim ^ Z t>L] = J_ . . 

L^ooQ[Z T ~Z t >L\Z t ] " T-t K ' 

9 Generalized inverse- Gaussian prior 

The generalized inverse-Gaussian (GIG) distribution is a three-parameter family of 



distributions on the positive half-line (see J0rgensen 23( or Eberlein & v. Hammerstein 



16J for further details). The density of the GIG distribution is 



f GIG (x; A, 5, 7) = 1 { .>o} (J) " 2K \ l5] xX ~ 1 ex P i-W 2 ^ 1 + 7 2 *)) • (60) 

Here K v [z] the modified Bessel function (see Abramowitz & Stegun |l|, 9.6]). The 
permitted parameter values are 

5>0, 7>0, if A > 0, (61) 

5 > 0, 7 > 0, if A = 0, (62) 

5 > 0, 7 > 0, if A < 0. (63) 

For A > 0, taking the limit 5 — > + yields the gamma distribution. For A < 0, 
taking the limit 7 — > + yields the reciprocal-gamma distribution — this includes the 
Levy distribution for A = —1/2 (recall that the Levy distribution is the increment 
distribution of stable-1/2 subordinators) . The case A = —1/2 and 7 > corresponds 
to the inverse- Gaussian (IG) distribution. If X has density ( 160]) then the moment 
/ifc = E[X fc ] is given by 

"' = W©' ^XeR,S>0,,>0, (64) 

r[A + *U2V k> _ x 



fi k =--{ r[A] V7 2 / and A >0, 5 = 0, 7>0, (65) 

00 k < —A 



H k = ^ r[-A] V 2 / and A < 0, 5 > 0, 7 = 0. (66) 

00 fe > —A 
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The following identity is useful (l|, 10.2.15]: 



K n+h [z} = J^-e- z J2(n + lj)(2zy\ forneN, (67) 

where (m, n) is Hankel's symbol, 

(m,n)= [ % J (68) 

n\ 1 [m + 2 — n\ 

The inverse- Gaussian (IG) process is a Levy process with increment density 

ct 1 / 7 2 / \ 2 \ 

9t (x) = 1 ^>0} ^^ a , 3 /2 eX P (^-^ ( X " P) J ■ ( 69 ) 

We see that qt(x) = J'gig(x] ~ §> d, l)- The fcth moment of qt(x) is 

™ {k) ~ ' ^ 7 e 7Ci f-V +2 Kk-iMjct], (70) 



V 7T \7 

for k > 0. Using (|67|) . the first four integer moments simplify to 

m?> = *, (71) 

7 

m| 2) = 4(l + 7ct), (72) 

m? } = 4(3 + 3 7 ct + 7 Vt 2 ), (73) 

(4) = C 5 + x + g 7 2 c 2 t 2 + 7 3 c 3 t 3N (74) 

r 

9.1 GIG terminal distribution 

We shall see that the GIG distributions constitute a natural class of a priori distribu- 
tions for ultimate loss. With 7 > and c > fixed, we examine some properties of a 
paid-claims process {£*t} with time-T density fciciz; A, cT, 7). The transition law is 

tT e dy I U = x] = J^fUy - x) dy, (75) 

^a \c 1 ^s{dy;x) 

TT^y\i sT = X } = - l - w - j , (76) 



where 



V>o(dz; £) = faia{z\ A, cT, 7) dz, (77) 



«pi-4 ^r-- 



^(d^; = (1 - t)1{ z >5} q _£/ z u/2 /g/o(«; a > cT > 7) dz. (7Z 
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Writing 



we have 



2-Y 



cTJ 2K X [ 1 VT}' 



(79) 



c 2 (T-t) 2 l 21 



( V(R:,) = M1-I)e-^- / z*+s " {z _ y)3/2 1, 



y 



. (T-ty i2- 



«(1 - £)e-^ / (* + y) A +^ 6 2 ' 2 dz 



T>^ / v^ ' »/ ~3/2 

./O 

«V^ 1 2 (T _ t) / A+ l 



r/ . yrfv-vP-QJ {z + y) x+ *q T _ t {z)&z. (80) 

Given £ t y = y, the best-estimate ultimate loss is 

U tT = i/j t (R;y) / z^ t (d^;i/) = -f— — - T — — -. (81) 

J y J (z + y) A+ 2q T „ t {z)dz 

9.2 Case A = -1/2 

When A = — 1/2 we have 

^(R;j/) 1 c(t - g) f 7 *((y- x )- c (t-s)/<yf 



= g t _.(y-x). (82) 

Hence {£tr} is an IG process. Note that in this case {£tr} has independent increments. 

9.3 Case A = n — \ 

We now consider the case where A = n — |, for n G N+. For convenience we write 

q i t k) (x) = f GIG (x;k-l/2,ct, 1 ), (83) 

Hence we have q\ (x) = qt(x). The transition density of {£tr} is then 
Mfr£ f (v x) - a (v . fo^z + vrgr-tWdz 

= q t -s(y-x) 1 - k = o[k) J~ l J . (84) 

it svf I _^ n / n x (n-fc) j. v ' 
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When n — 1 this is 



-jt-s{y-x) = q t -s{y-x)- 






c(t-s) ^ im 



7X + c(T — s) 



c(t-s) \ (h 



Thus the increment density is a weighted sum of GIG densities. We shall now derive 
a weighted sum representation for general n. We can write 

/■oo 

(z + y) n q T -t(z) dz = / ((z + x) + (y-x)) n q T -t(z)dz 
o Jo 

= E (?) (y - x ) n ~ k jf(* + x ) k «r-*w d ^ 

^(ij^-^tf^-^ (86) 

fc=o ^ ' i=o ^' 

Then we have 



M^y), ,.. ^_„ ,.. ^JT^ + ^-^dz 



-ft- s {y-x) = q t - s (y-x)- 



- k ^ k ^mL fc_ fV' 



/ _ v ^Lo 0(2/ ~ x )"~ fc Ej=o (j) '"■/■-/ 
- q t - 8 {y x) 

E fc =o U m T- '■' , 
However, when k G No, 

z k q t -s(z) z k f GIG (z; -1/2, c(t - s), 7) 



& 



(fc) 



(z) " / Cffl (z;fc-l/2,c(i-s),7) 

c(t - s)\ h K k ^y 2 [yc(t - s)] 



7 / K l/2 [jc(t - s)} 

(fc) 



Hence we have 

\n—k „ /„, „\ _ ^,(™- fc ) „(. n ~ k ) 



(y - x) n ~ h q t _ s (y - x) = m\T s K) C % - x). (89) 

Using the identity (|89|) . ( 157)1 can be expanded to obtain 

^ s(i ' ' fc=0 



where 



,( fc ), 



,.(«-*) v* . r fc w( fc 7V 



G)m£* ; EL, ( •) 



l T-t 



EL,CW-7^- 



(91) 



J 3=V \J 



Notice that w st (x) is a rational function whose denominator is a polynomial of order n, 
and whose numerator is a polynomial of order k < n. Thus the transition probabilities 
°f {£tr} depend on the first n integer powers of the current value. The conditional law 
of the ultimate loss is 



n„(0) 



^(dy;&r) y n QT-s(v ~ 6r 



^(R;W EJU&"fc 



(ra-fc) 
s 



dj/. 



(92) 



Ml (fc). 



We can verify that E^o^st V^) = ^ using the fact that IG densities are closed 
under convolution. We have 



Qt-s{z) 



QT-t{y)Qt-s(z -y)dy, for < s < t < T. 



For fixed n G N+, we then have 

n 

E 



k=0 



n \ (n-fc) fc 



(z + a;) n g T _ s (2;)d2; 



{z + x) n J q T ^ t {y)qt- s {z - y) dy dz 
o 

"OO 

q T -t(y) I (z + x) n q t _ s (z - y) dz dy 



Qr-t(y) / (z + y + %Y Qt-s(z) dz dy 
o Jo 



Qt-a(y) 



k=0 



n 



J2\ >. )'">-» 



(n—k) 



(y + x) 



dy 



E 

k=0 



n \ (n~k) 



(y + x) k q t - s (y) dy 



eO^eQ 

k=o v 7 j=0 VJ/ 



m T _ t x 



which gives 



E 

k=0 



(>fc fc ^U-K-7^ 



,(n-J) 



e; =0 (>^' 



i. 



(93) 



(94) 



(95) 
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9.4 Moments of the paid-claims process 

The best-estimate ultimate loss simplifies to 

+ 1 (n+l\^(n+l-k) c k 



U = Z^fc = I k ) " l T-t UT /ggX 

2^k=0 \k) m T-t str 



For example, when n = 1 we obtain 



c(T - t)(l + 7 c(T-t)) + 2 7 2 c(^- t)itT + 7 3 e t 2 



tt H- 1 _ L )\^ -r- 7H- 1 ^ ^ -*7 H- 1 ^«TT7^ ff , , 

^*r 7777^ — TT-TT^T ■ I 97 J 



By similar calculations, we have 

v-^'i+m (n+m\ (n+m-k) j-k 

Wfifr | 6t] = ^ fc= ° l fc j r "l. formed, (9? 

Z^fc=0 \k)" L T-t StT 

and 

(n-k) c k 



E 



g2 a2 fTT 



;*r 



Z^fc=o \k) m T-t str 



for < a < 7, where 7 = a/7 2 — a 2 , and mj is the /cth moment of the IG distribution 
with parameters S — ct and 7 = 7. 

10 Exposure adjustment 

We have seen that 

E[6t] = |e[C/ t ]; (100) 

thus in the model so far the development of the paid-claims process is expected to be 
linear. This is not always the case in practice. In some cases the marginal exposure is 
(strictly) decreasing as the development approaches run-off. This manifests itself as 

q2 

o^nctT] < 0, (101) 

for t close to T. A straightforward method to adjust the development pattern is through 
a time change. We describe the marginal exposure of the insurer through time by a 
deterministic function e : [0,T] — > R+. The total exposure of the insurer is 

T 

e(s)ds. (102) 
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We define the increasing function r(t) by 

. , L e(s) ds 
r(t) = T^A^ . (103) 

Jo £ ( s ) ds 

By construction r(0) = and r(T) = T. 

Now let r(t) determine the operational time in the model. We define the time- 
changed paid-claims process {Q T } by 

e[ T = e(r(t),T), (104) 

and set the reserving filtration to be the natural filtration of {Q T }- Then we have 

rt i \ j 

E[&] = J ° 7; * Ep7 T ] (105) 

Jo e(s)ds 



and 



ot J o £ (s)ds 



10.1 Craighead curve 



Craighead [15J proposed fitting a Weibull distribution function to the development 
attern of paid claims for forecasting the ultimate loss (see also Benjamin & Eagles 
. In actuarial work, the Weibull distribution function is sometimes referred to as the 
'Craighead curve'. To achieve a similar development pattern we can use the Weibull 
density as the marginal exposure: 

e(t) = -(x/a) b - 1 e- {x/a) \ (107) 

for a, b > 0. Then the time change r(t) is the renormalised, truncated Weibull distri- 
bution function 

1 _ e -(t/a) b 

r ® = T 1 _ e -(W ( 108 ) 

See Figure [2] for plots of this function. When b < 1, r'(t) is decreasing. Under such a 

time change, the marginal exposure is decreasing for all t e [0,T]. When b > 1, r'(£) 

achieves its maximum at 

* fb-l\ 1/b , , 

t* = a I — ) , (109) 

and r'(t) is decreasing for t > t*. Thus, if T > t* then the marginal exposure is 
decreasing for t E [t*,T]. If T < t* then the marginal exposure is increasing for 

te[o,T]. 
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Figure 2: Plots of the truncated Weibull time change for various parameters, and with 
T = 1. The expected paid-claims development of the model will have the same profile 
as r(t) (scaled by E[C/r])- Hence, under one of the above time changes, when t is close 
to T the marginal exposure falls (i.e. J^-E[£ 4 T T ] < 0). 

11 Simulation 

We consider the simulation of sample paths of a stable-1/2 random bridge. First, we 
can generalise (1H]1 to 



[^ s ^,T)\^8,T)=y^(t,T) = z] 



law 



1/ x / z 

y + -(z - y) H ; 

2 y> \ v / c 2 (t _ a )2 /(jz _ y) + Z 2 



(110) 



where < s < t < T, and Z is a standard normal random variable. We can then 
generate a discretised sample path {£(£j,T)}^™ , where tj = iT2~ n , by the following 
recursive algorithm: 

1. Generate the variate £(T,T) with law is, and set £(0, T) = 0. 

2. Generate |(f ,T) from |(0,T) and C(^,^) using identity ffTTU]) . 

3. Generate ^(?,T) from ^(0,T) and ^(^,T), and then generate £(*f,T) from 
|(f, T) and |(T,T). 

4. Generate |(f ,T), |(f ,T), |(f ,T), |(f ,T). 

5. Etc. 

See Figure [3] for example simulations. 
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Figure 3: Simulations of the paid-claims process {£,tr} (bottom line) and the best- 
estimate process {Utr} (top line). Various values of the activity parameter c are used. 
A priori, the ultimate loss Ut has a generalized Pareto distribution (GPD) with density 
Igpd(x) = l{x>i} (l + ^-) ■ (This is the GPD with scale parameter a — 1, location 
parameter \x — 1, and shape parameter £ = 1/4.) 

12 Multiple lines of business 

We shall generalise the paid-claims model to achieve two goals. The first is to allow 
more than one paid-claims process, and allow dependence between the processes. The 
second is to keep the dimensionality of the calculations low with a view to practicality. 
The following results can be applied to the modelling of multiple lines of business or 
multiple origin years when there is dependence between loss processes. 

12.1 Two paid-claims processes 

We consider a case with two paid-claims processes, but the results can be extended 
to higher dimensions. In what follows, we set ft(x) = ft(x) as given by (jl]), and 
f£ T {x) = ftT^x) as given by (fTTj) . Here we have introduced the superscript to emphasise 
the dependence on c. Let {S(t, T*)} be a stable-1/2 random bridge with terminal 
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density p(z) = u(dz)/dz, and with activity parameter c. Fix a time T < T*, then 
define two paid-claims processes by 

$ = S(t,T*) (0<t<T), (111) 

£§> = k 2 S{Xt + T, T*) - k 2 S{T, T*) (0<t<T), (112) 

where A = T* jT — 1, and k = c^/(cX) for some c 2 > 0. The density of £^f is given by 

^\x) = r T {x)r f ^} z - x) P {z)dz 

Jo Jt*\ z ) 

oo 



f^(x;z)p(z)dz, (113) 

and the density of £^ is 



'o 

.(2) 



p^\x) = k- 2 f^_ T (k- 2 x) r m * ;\ 2x) p{z) dz 

Jo Jt*\ z ) 

= k^f^_ T (k- 2 x) r mk ~**~ k ~ 2x) p(k-*z) dz (114) 

Jo Jt*\ K Z > 

= k- 4 f^_ T)T ,{k- 2 x]k- 2 z)p{k- 2 z)dz (115) 

Jo 

= k~ 2 f^\^ T ,{x-z)p(k- 2 z)dz. (116) 

Jo 

Here (11141) follows after a change of variable, (I115p follows from the definition of frr(y, z) 
given in ffTUj) . and (I116p follows from the functional form of frr(y,z) given in ffTB . 
It follows from the dynamic consistency property that {Qtj * s a stable- 1/2 random 
bridge with terminal density p^(z) and activity parameter c. Using the dynamic 
consistency property and the scaling property of stable-1/2 bridges, one can show that 
{Q T } is a stable-1/2 bridge with terminal density p^(z) and activity parameter c 2 . 
The conditional, joint density of (Q T , , k~ 2 Q T ') is 



r(i) c J., l.-2c(2) 



r(i) - ~ L-2 t (2) 



^ J G djfc, k- 2 e t r' G dy 2 £# = si, ^£r = *a 

/T>-(l+A) f ( 2: -(^+l/2)) , . , 1 „ , . , „ , , , 

— — - — rrp{z) dz )■ f t _ s {yi - xi) dyi / A(t _ s) (y 2 - x 2 ) dy 2 , 

ci+x 2 JT*-(l+X) s y Z \ x l-r x 2)) J 

(117) 
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for < s < t < T. Then we have 



r(i) 



2t(2) 



Qt' + k-Ttr' G dy 



r(i) 



2 t (2) 



SsT ~~ X l' ^ SsT — x 2 



./• 



i+x 2 JT*-(1+X) 



T * ^ (t t ; ^ p(^) dz f f(i+x)(t- S )(v - fa + a*)) dy 



2=xi+2;2 

t-(l) 



/(i+A)(t- s ),T*-(i+A) s (> - fa + ^2); V ~ fa + x 2 ))p(z) dz \ dy; (118) 



-2t(2) 



(i) 



and, given ^ T ; = x\ and fc ^ T = x 2 , the marginal density of Q T ' is 



l/i >->■ 



/ t C _s,T*-(l+A) S (l/l -Xl',Z- fa + X 2 ))pfa dz, 



z=x\+x<z 



119) 



and the marginal density of fc 2 £ 4T is 



'A> ! ?■ / /a(<- s ),t*-(i+a) s (2/2 -x 2 ;z- (zi + x 2 ))pfa dz. (120) 



12.2 Correlation 

The a priori correlation between the terminal values is well defined when the second 
moment of v is finite. The correlation can be used as a tool in the calibration of the 
model. Assuming that W,[S(T*,T*) 2 } < 00, the correlation is defined as 



E 



t(i) <r(2) 



-E 






E 



t(2) 



E 



^vttJ 


-E 




■)(. 


1 yrr) 


-E 


>(2)" 
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We shall calculate each of the components of (11211) separately. First, we have 



Noting that 



E 



t(i) 



E[£(T,T*)] = ^E[S(T*,T*)]. 



,(2) 



^ = k\S(T\T*) - S(T,T*)) 
l = k 2 S(T*-T,T*), 



122) 



(123) 



we have 



E 



t(2) 



fc 2 E[S(T*-T,T*)] 

p(\--^W(r*,r*)]. 



(124) 
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f(l) 



r(2) 



The second moments of £yj, and ^^ follow from (H3]) , and are given by 



E 



(e 



(i) 

TT 



T 



E [S(T*,T*) 2 ] - (T* -T)C T *, 



(125) 



and 



E 



(6 



(2) 



T 



T» 



fc 4 1 E [S{T\T*f] - k A TC T * 



where 



c i T * 



C T , = cV^ / z 3/2 e^~ $ [-cT*z~ 1/2 ] p{z) dz. 



(126) 



(127) 



The final term required for working out the correlation is the cross moment. This is 



E 



t(l) A 2 )' 

STT ?TT 



A; 2 E [S(T, T*) (S(T*, T*) - S(T, T*))] 



= k 2 E [S(T, T*)S(T*, T*)] - k 2 E [S(T, T*) 2 ] . 
The first term on the right of (I128p is 



128) 



oo /-oo 



A; 2 E[^(T,T*),S(T*,T*)] = fc- / / xy 

'o </o 

*7 / 



#.(v) 

xyfT T *(x;y)dxp(y)dy 



dxp(y) dy 



rp /*0O 

k 2 — / y 2 p(y)dy 
k 2 ^-E[S(T*,T*) 2 ]. 



The second term on the right of fll28[) is given by (11251) . Hence we have 



E 



r(i) t( 2 ) 



^^ =A; 2 (T*-T)C T ,. 



(129) 



(130) 



The expression for the correlation follows by putting equations (I122p . (I124p . ( I125p . 
( TT26]) . and ffT30l together. 



12.3 Ultimate loss estimation 

We now estimate the terminal values of the paid-claims processes. At time t < T, the 
best-estimate ultimate loss of {Q T } (or, indeed, {Q T }) depends on the two values £ tT 
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,(2) 



and £ tT . The best-estimate ultimate loss of {Q T } is 



(i)i 



£/, 



(i) 



E 






t(l) 



J/1, C, t y 



:r 2 



E [S(T, T*) | 5(t, T*) = xi, S(T + At, T*) - S(T, T*) = k~ 2 x 2 ] 

E [S(T + At, T*) | S(t, T*) = xi, S{T + At, T*) - 3{T, T*) = k~ 2 x 2 ] - k~ 2 x 2 

E [S(T + At, T*) | 5(t, T*) = xi, 5((1 + A)t, T*) - S(t, T*) = k~ 2 x 2 ] - k~ 2 x 2 

(131) 

E [S(T + At, T*) | 5((1 + A)t, T*) = x x + k~ 2 x 2 ] - k~ 2 x 2 (132) 

' "" ' (E [S(T*,T*) | S((l + A)t,T*) =xi + k~ 2 x 2 ] - k' 2 x 2 ) 



m* 



;i + A)t 

T* - (T - t) 

+ T*-(l + A)t Xl ' 



(133) 



Equation f)13ip holds since reordering the increments of an LRB yields an LRB with 
same law, (J132J) follows from the Markov property of LRBs, and f)133p follows from 
El). We also have 



POD 

E [S(T*, T*) | 5((1 + A)t, T*) = xi + A- 2 :r 2 ] = / zp 4 (z) dz, 

JO 



(134) 



where 

p t (z) = l{ 2 >xi+fc- 2 x 2 }-^ _1 



fxi + k 2 x 2 ) 



3/2 



c 2 / (T* - (1 + A)t) 2 T* 2 \ . 



and fC is a constant chosen to normalise the density. Similarly, the best-estimate 
ultimate loss of {Q T } is 



U< 



(2) 



fT 



fc 2 T* (T .? (E [5(T*,T*) | 5((1 + A)t,T*) = xi + k~ 2 x 2 ] - Xl ) 



rrjf 



X)t 



+ 



T-t 



T*-(1 + A)t 



x 2 . (136) 



To compute both UJ: T and £/ tT we need to perform at most two one- dimensional in- 
tegrals (the integral we need is (11341) . but we note that pt(x) includes a normalising 
constant K — which is found be evaluating a second integral). We are saved the com- 
plication of performing double integrals. 

To extend these results to higher dimensions we can split the 'master' process 
{Str} into more than two subprocesses. Regardless of the number of subprocesses 
(i.e. paid-claims processes), all of the best-estimate ultimate losses can be computed 
by performing at most two one-dimensional integrals. This makes such a multivariate 
model highly computationally efficient. 
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